Hybrid improper ferroelectricity: a mechanism for controllable 
polarization-magnetization coupling 
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First-principles calculations are presented for the layered perovskite Ca3Mn207. The results 
reveal a rich set of coupled structural, magnetic and polar domains in which oxygen octahedron 
rotations induce ferroelectricity, magnetoelectricity and weak-ferromagnetism. The key point is 
that the rotation distortion is a combination of two non-polar modes with different symmetries. We 
use the term "hybrid" improper ferroelectricity to describe this phenomenon and discuss how control 
over magnetism is achieved through these functional antiferrodistortive octahedron rotations. 
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The utility of multiferroics for low-power electronic de- 
vices stems from the possibility for electric-field control 
of magnetism at room temperature ^0]. A challenge 
that has so far not been overcome is to identify a stable, 
single phase multiferroic material in which the magneti- 
zation can be deterministically switched 180°. A large 
electrical polarization strongly coupled to the magneti- 
zation is generally thought to be a key requirement 0]. 

Magnetically-driven improper ferroelectrics, such as 
TbMnOa, are materials in which a spontaneous polar- 
ization arises due to symmetry-breaking by a spin in- 
stability [3 |6j . These materials naturally have a strong 
coupling between magnetism and the polarization, but 
the polarization is too small for device applications. In 
known multiferroic materials with a large electrical polar- 
ization, such as BiFeOa f?! , the ferroelectricity is proper, 
originating from a zone-center polar lattice instability, 
as in the prototypical perovskite ferroelectric PbTiOs. 
However, except in a few special cases that satisfy restric- 
tive symmetry criteria, the polar instability in a proper 
ferroelectric does not break the right symmetries to turn 
on a nonzero magnetization and therefore does not satisfy 
the criteria of the only known mechanism that enables 
the electric field switching of the magnetization (S] . 

It is desirable to identify a more general mechanism 
- applicable to a large class of materials, for example, 
the ABO3 perovskites - whereby ferroelectricity and fer- 
romagnetism are induced by the same lattice instability. 
Octahedron rotations, ubiquitous in perovskites and re- 
lated materials, are natural candidates for this lattice 
instability as they are known to strongly couple to mag- 
netic properties [9HTT]. Unfortunately, such distortions 
(or combinations of distortions) in simple perovskites are 
not polar and therefore do not induce ferroelectricity. Re- 
cently, however, Bousquet, et al. [T^] made a key discov- 
ery that by layering perovskites in an artificial superlat- 
tice, e.g., (SrTi03)/(PbTi03), a polarization can arise 
from the coupling of two rotational modes. Taking ad- 
vantage of this mechanism to realize a strongly coupled 
multiferroic is a challenge. 

In this Letter, we demonstrate how octahedron ro- 
tations simultaneously induce ferroelectricity, magneto- 
electricity, and weak ferromagnetism in a class of natu- 



rally occurring (AB03)2(AO) layered perovskites. The 
key point is that the polarization, P, arises from a ro- 
tation pattern that is a combination of two non-polar 
lattice modes with different symmetries, P '^ TZiTZ2, 
as in Ref. [T^ but here rotations TZi and TZ2 addition- 
ally induce magnetoelectricity and weak ferromagnetism 
respectively. We use the term "hybrid" improper fer- 
roelectricity to describe this ferroelectric mechanism (in 
loose analogy to improper ferroelectricity [13] such as in 
YMn03 [14]) in order to generalize the idea to include 
cases where the two distortion patterns do not necessar- 
ily condense at the same temperature. This mechanism 
has no impediment to room temperature operation and 
in fact opens up entirely new classes of materials in which 
to search for strongly-coupled multiferroics. Our results 
show a rich set of coupled structural, magnetic and polar 
domains and suggest the possibility to switch between 
magnetic domains with an electric field. 

We have identified (CaB03)2CaO, with B=Ti ^E\, 
Mn [ini [17], as two materials that display hybrid im- 
proper ferroelectricity. It is significant that they occur in 
nature in bulk, forming in the Ruddlcsdcn-Popper (RP) 
homologous family with general formula A„+iB„03„+i. 
Any given member of the RP series consists of ABO3 
perovskite blocks stacked along the [001] direction with 
an extra AO sheet inserted every n perovskite unit cells. 
For Ca3Mn207 (n = 2) the experimental picture of 
the sequence of phase transitions from the paraelec- 
tric lA/mmm phase to the ferroelectric A2iam phase 
is not clear. Two possibilities have been proposed: (1) 
lA/mmm — >■ Cmcm — >■ A2iam, and (2) a direct transi- 
tion from lA/mmm — >■ A2iam. Additionally, it has been 
shown to display weak ferromagnetism |17l I18| . To our 
knowledge, Ca3Ti207 has only been reported in the po- 
lar A2iam structure. In the remainder of this Letter, we 
focus mainly on the magnetic compound, Ca3Mn2 07, as 
a prototype of this class of strongly coupled multiferroics. 

First-principles calculations were performed using 
DFT using PAW potentials within LSDA+U [19 as im- 
plemented in VASP [^Ul - ES] . All calculations were re- 
peated with the PBEsol functional, which provides an 
improved description of structural parameters; there was 
no qualitative change in any of our results. We used U = 





(a) 







(b) 




"%. / 




0.05 


\ / 


0.0015 


-0.1 


Vy 


> 

< 


0.15 


4 ■ \/ 


0.0005 


-0? 


X3,« 







0.2 0.4 0.6 0.8 1 1.2 1.4 
Q(A) 



0.06 -0.03 0.03 0.06 
Q(A) 



FIG. 1: Ca3Mn207 structure and rotation distortions, a, 
The A2iam ferroelectric ground state structure. Blue (large) 
spheres correspond to Ca ions, b Schematic of the atomic 
displacements corresponding to the X^ rotation. The dashed 
square denotes the unit cell of the M/mmm parent struc- 
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FIG. 2: Energy surface about paraelectric M/mmm in 
Ca3Mn207. Change in energy per formula unit as a function 



ture. c, Schematic of the X^ tilt mode. All axes refer to the of the amplitude of a, the XJ tilt and X3 rotation modes and 



coordinate system of the H/mmm, parent structure. 



4.5 eV and Jh=1 eV for the Mn-ion on-site Coulomb and 
exchange parameters respectively. Where noted, non- 
collinear calculations with L-S coupling were performed. 
We used a 600eV plane wave cutoff, a 4x4x2 Monkhorst- 
Pack mesh, which we checked for sufficient accuracy. 

In the polar A2iam structure, the oxygen octahedra 
are significantly rotated and tilted with respect to the 
lA/mmm structure, as shown in Figure nl The polariza- 
tion, in the a;j/— plane by symmetry, is found from our 
first-principles calculations to be large, P « 5 /iC/cm^ 
{P « 20 AiC/cm2 foj. CaaTiaOy). Group theoretical 
methods show that A2iam is related to lA/mmm by 
three distinct atomic distortions: a polar zone-center 
mode transforming like the irreducible representation (ir- 
rep) T^ , and two zone-boundary modes at the X (1/2, 
1/2, 0) point - an oxygen octahedron rotation mode with 
irrep X2 and an oxygen octahedron tilt mode with irrep 
X^ . Note that X2 ® X^ establishes the A2iam space 
group, a zone-center polar instability is not required. 
Hence, it is possible to reach the ferroelectric state by 
means of a combination of rotations and tilts only. 



b, the polar Fg distortion, c and d. Polarization in the pres- 



ence of a hybrid order parameter, Qx2a- 
the differences in scales between panels. 



=AQ„-Q„+. Note 



resulting in a ground state with four structural domains. 
Figure [2b, d shows how the polarization arises from 
a coupling to a "hybrid" order parameter QX23 = 
Qx+Qx^ ■ ^^ ^^^ absence of rotation and tilt distor- 
tions ((3x23 =0)1 the polarization has a single minimum 
at P =0. As Qx23 increases, the polarization never be- 
comes unstable. Rather, the minimum shifts to a non- 
zero value. The result of increasing QX23 is analogous to 
the effect of turning on a finite electric field, just like in 
the classic case of improper ferroelectricity |13j . Further- 
more, when Qx23 ¥" Oj the polarization is linear about 
zero (Figure^), a direct indication of improper coupling 



j^ = aPQx+Q 



We project out the contribution of each X^ , X2 , and 
r^ mode to the A2iam ground state structure and calcu- 
late from first principles the T = energy surface around 
the lA/mmm reference structure. Figure^ shows the to- 
tal energy as a function of the amplitude of the distortion 
for the individual rotation ((5^+), tilt {Qx-)^ 8'iid polar 
(Qp-) distortions. Relatively large energy gains can be 
seen within a characteristic double-well potential for the 
rotation and tilt distortions, whereas the polar contribu- 
tion is stable, as shown in Figure[2]3. Additionally, Figure 
3 shows that the combination of Qx+ plus Q x- lowers 
the energy - even in the absence of the Qp- distortion ~ 



'^x: 



between the polarization, rotations, and tilts. 

These observations suggest that the rotation mode and 
the tilt mode are the primary modes driving the transi- 
tion to the ferroelectric A2iam phase. This single distor- 
tion pattern, QX231 is the hybrid improper mode. In con- 
trast to proper and conventional improper ferroelectrics, 
more than one lattice distortion may switch the polar- 
ization in a hybrid improper ferroelectric. Symmetry im- 
plies, and our calculations confirm, that the polarization 
reverses by either switching a X^ distortion or a X^ 
distortion but not both, resulting in two polar domains. 

As is often the case in perovskites and related mate- 
rials, octahedron rotations directly couple to the mag- 
netic ordering. We determined from first principles the 
magnetic ground state of A2iam Ca3Mn207 with po- 
larization along [010] to be antiferromagnetic (G-type 
within the perovskite bilayer). The spins point along 
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FIG. 3; Four types of structural domains and possible paths 
for electric-field switching of the magnetization in Ca3Mn207. 
a, Unstrained (bulk): Starting from a fixed X^ tilt domain, 
switch polarization by switching X2 rotation: Pq — >■ — Pq and 
Mo — ^ Mq. b, Compressive biaxial strain: Start from a fixed 
X2 rotation domain, switch polarization by switching X^ 
tilts: Po — >■ — Po and Mn — )■ — Mq. The energy change is per 
formula unit. 



[001] due to crystalline anisotropy. Additional spin-orbit 
interactions give rise to a net spin-canted moment of 
M ss 0.18/z_B per unit cell (4 spins) along [100]. These 
results are consistent with previous experiments. Note 
that the magnetic point group, 2'mm' , allows for a lin- 
ear magnetoelectric effect [21] , which symmetry indicates 
is induced by the X^ rotation distortion. Application of 
Dzyaloshinskii's criteria [25l [26] shows the canted mo- 
ment is the result of the X^ tilt distortion. Indeed, if we 
compute the magnetic ordering in the Cmcm (Cmca) 
structure obtained by freezing in the X^ ( X2) mode 
alone, we find M w 0.22/iB (M = O/is) per unit cell, 
which reverses with reversal of this octahedral tilt ^7\ 
resulting in two magnetic domains. 

So which distortion, the rotation or the tilt, reverses in 
an electric-field switching experiment? Although polar- 
ization switching in a ferroelectric is a complex, dynamic 
process, we may gain some insight by examining the in- 
trinsic energy barriers between domains. As shown in 
Figure |3^ the lowest energy pathway to switch the direc- 
tion of the polarization is along the X2 switching path. 
In this process the magnetization does not reverse its 
sign. The linear magnetoelectric effect, however, is in- 
duced by the X2 distortion as mentioned. This electric- 
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FIG. 4: Efi'ect of epitaxial strain on the unstable zone- 
boundary modes in Ca3Mn207. Change in energy per for- 
mula unit w.r.t. paraelectric M/mmm, as a function of the 
amplitudes of the X^ and X^ modes under a, 1.5% tensile 
and b, 1.5% compressive biaxial strain. The lines are a guide 
to the eye (similar results were also found for Ca3Ti207). 



field tunable oxygen rotation distortion may lead to an 
enhanced magnetoelectric effect. Such calculations are 
beyond the scope of this Letter, but future theoretical 
and experimental studies should make this clear. 

Oxygen rotations in perovskites are known to respond 
strongly to pressure and epitaxial strain. Figure[3b shows 
the energy landscape around the lA/mmm paraelectric 
structure at 1.5% compressive strain. Now the lowest en- 
ergy pathway to switch the polarization to a symmetry 
equivalent state is along the X^ switching path, which as 
previously discussed, switches the direction of the spin- 
canted moment. Therefore, for an epitaxial thin film 
compressively strained in the A2iam phase, we predict 
that switching the direction of the polarization with an 
electric field will switch the direction of the equilibrium 
magnetization by 180°. As in all experiments to date 
based on the linear magnetoelectric effect, however, a 
single antiferromagnetic domain must be annealed and 
maintained throughout the experiment. 

This observation of tuning the intrinsic energy barriers 
between the domains with strain can be understood from 
well-known simple physical considerations |2H]- Figure [4] 
shows the behavior of the rotation and tilt distortions in 
Ca3Mn207 under 1.5% biaxial tensile and compressive 
strains. Figure [4k shows that the energy lowering of the 
X2 rotation is strongly reduced under tensile strain com- 
pared with the unstrained state shown in Figure [2k. Un- 
der compressive strain, the opposite behavior occurs: the 
X2 mode is strongly favored, lowering the energy even 
more than the X^ tilt mode, as shown in Figure kp. We 
emphasize that strain does not induce, enhance or medi- 
ate the polarization-magnetization coupling responsible 
for the electric-field switching of the magnetization in 
these materials. In the specific case of Ca3Mn207, strain 
merely alters the energy landscape around the lA/mmm 
paraelectric structure, biasing one hybrid improper fer- 
roelectric switching pathway (the one that simultane- 
ously switches the direction of the magnetization) over 



the other. However, there is no fundamental reason why 
electric-field switching of the magnetization cannot be 
observed in bulk hybrid improper ferroelectrics. 

In addition to the measurements already suggested, 
spatially resolving the structural, polar, and magnetic 
domains, e.g., optically [22], should prove the coupHng 
physics discussed even in bulk Ca3Mn207. It would also 
be of interest to understand the phase transition sequence 
from the high- symmetry paraelectric lA/mmm phase to 
the low-symmetry ferroelectric A2iam phase, which is 
impossible to prove or disprove from T=0 calculations. 
In bulk, two possibilities have been proposed as we pre- 
viously discussed [SD]. Path (1), having an intermediate 
Cmcm phase, is consistent with Landau theory and with 
our calculated hierarchy of structural distortions as dis- 
played in Figure [3^. If Path (2) turns out to be correct, 
however, this indicates that the two distortions making 
up the hybrid order parameter condense at the same 
temperature, exactly analogous to the recent discovery 
of Ref. [121 It is unlikely to be a coincidence that both 
PbTiOa/SrTiOa [H] and (CaMn03)2/CaO are layered 
perovskites. Regardless of the actual path taken - or 
even if a paraelectric structure is realizable in the phase 
diagram of the system - our conclusions on the coupling 
of rotation/tilt distortions to ferroelectricity and mag- 
netism remain unchanged. 

Note that the temperature scale of the hybrid improper 
ferroelectric mechanism is set by the structural distor- 
tions, which commonly occur above room temperature 
(~500-600K in CaaMuaOr 1310- The 'limiting' temper- 
ature in the case of Ca3Mn207 is the Neel temperature, 
T/v ^115K [171 [ISj and as such this prototype material is 
not ideal. There is, however, no fundamental reason why 
one couldn't discover (or design) a hybrid improper ferro- 
electric with a Neel temperature above room temperature 
(ongoing investigations into the connection between lay- 
ering and hybrid improper ferroelectricity have realized 
initial materials design rules [32] ) , so the mechanism has 
no impediment to room-temperature operation. This ap- 
proach therefore shifts the challenge of discovering a mul- 
tiferroic in which the magnetization can be controlled by 
the electrical polarization to the more familiar problem 
of designing a room temperature antiferromagnet. 

In summary, we have introduced the term hybrid 
improper ferroelectricity to describe a state in which the 
polarization is induced by a complex distortion pattern 
consisting of more than one octahedron rotation mode. 
We have shown how this mechanism opens up new 
avenues to pursue strong polarization-magnetization 
coupling, as alluded to in Ref. [12] and realized in this 
Letter, and discovered a new class of materials. We 
hope our Letter inspires further work in this intriguing 
field of materials with functional oxygen rotations and 
tilt distortions. 
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AUXILIARY MATERIAL 

Landau theory to describe the polarization 
magnetization coupling depicted in Figure 3. 

Note that in the main text we use the tetragonal, 
lA/mmm structure (which has one formula unit per cell) 
as the paraelectric reference structure to discuss the hy- 
brid order parameter. To make clearer the polarization- 
magnetization coupling, in this auxiliary material we use 
two different orthorhombic reference structures (which 
each have two formula units per cell) corresponding to 
the two different possible saddle points shown in Figure 3 
of the main text. Therefore due to the subsequent zone- 
folding, the atomic distortions that make up the order 
parameter that we refer to below as the polarization, V, 
differ from those that make up the polar order parameter, 
P, from lA/mmm defined in the main text. This does not 
change the symmetry of the coupling of magnetization to 
the physical polarization, and hence does not change the 
conclusions that octahedron rotations induce ferroelec- 
tricity, magnetoelectricity and weak- ferromagnet ism. 

In the antiferromagnetically ordered phase, we perform 
an energy expansion in the magnetization and in an order 
parameter, P, that has similar transformational proper- 
ties as the polarization about the saddle point along the 
X^ switching path, i.e., at Qx- 7^ & Qx+ — 0- This 
corresponds to an orthorhombic Cmcm paraelectric ref- 
erence structure, for which symmetry analysis and Fig- 
ures [5^ and b show the ordered antiferromagnetic vector 
'Lq — Si ~ S2 — S3 + Si being even with respect to space 
inversion, i.e., ILo = 54 — 6*3 — iS'2 + 5'i — -t-Lg. Note that 
we confirmed that Cmcm is the saddle point structure by 
performing full structural relaxations from first princi- 
ples within this Cmcm tilt only structure and within the 



where 7 is a coupling coefficient. A standard Landau 
approach [1] leads to 



Cmca rotation only {Qx- = & Qx+ ¥" 0) structure. 
The tilt only structure is indeed closer in energy to the 
A2iam ground state. Here, in this antiferromagnetic- 
paraelectric structure symmetry allows the magnetiza- 
tion, M, to couple to even powers of V, where the lowest 
order coupling is of the form 



Mo 

Am 



7^0 + 0{VS) 
j2{Lo'Po)A-p, 



(1) 
(2) 



where A-p and Aj\/ are small changes about the equilib- 
rium values of the magnetization (Afg) and Vq due to 
an applied electric field. This analysis makes clear that 
switching the physical polarization with an electric-field 
to its symmetry-equivalent state by following the hybrid 
order parameter along the X2 switching path does not 
switch the direction of the magnetization. The X2 rota- 
tion mode, however, does induce the linear magnetoelec- 
tricity as seen in Eq.2. 

Similar to the previous analysis, we perform an en- 
ergy expansion in V and magnetization about the saddle 
point along the X^ switching path, i.e., Qx- = & 
Qx+ 7^ 0. This corresponds to an orthorhombic Cmca 
paraelectric reference structure where symmetry analysis 
and Figures [5^, c, and d show the ordered antiferromag- 
netic vector Lq ~ Si ^ S2 — S^ + 6*4 is now odd with 
respect to space inversion, i.e., ILq = ^3 — ^4 — 5i + ^2 
= -Lq. We again confirmed that Cmca is the saddle point 
structure by performing full structural relaxations. Here, 
in this antiferromagnctic-paraelectric structure symme- 
try allows M to couple to odd powers of V, where the 
lowest order coupling is of the form 



Ei. 



^LoVM, 



where ^ is a coupling coefficient. Again, using a standard 
Landau approach [1] leads to 



Mo ~ ^LoVo 
Am ~ ^LqA-p, 



(3) 
(4) 



where we now see that the equilibrium magnetization is 
proportional to the equilibrium polarization [2]. It van- 
ishes at the paraelectric saddle point where Vq — 0, and 
is directly proportional to the strength of the ferroelec- 
tric distortion, i.e., X^ induces weak-ferromagnetism. 
Therefore switching the physical polarization with an 
electric-field to its symmetry-equivalent state by follow- 
ing the hybrid order parameter along the X^ switching 
path does indeed switch the direction of the magnetiza- 
tion. 
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FIG. 5: Effect of ^^3^ and X2 on antiferromagnetic vec- 
tor within a bilayer of Ca3Mn2 07. Schematic showing 
inversion centers, x, with respect to the spins and antiferro- 
magnetic vector liQ = Si — S2 — S3 + S4 for a and c, undis- 
torted octahedra {M/mmm) and b, X^ {Cmcm) and d, X2 
(Cmca) distorted octahedra. Note in the bottom figure that 
the inversion center along [001] between, e.g., Si and S3, is 
preserved. All numbers refer to the unique spins located in 
the center of the octahedra. 



